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INTRODUCTION 



In the past few years there has been significant progress in the study of the chiral- 
continuum extrapolation problem for lattice QCD with different types of fermions (see, for 
example, Refs. [H 0,0,0]). Extrapolation of lattice QCD simulation results to the physical 
(light) quark masses is a nontrivial task due to the nonanalytic variation of hadronic prop- 
erties with quark masses which is a consequence of spontaneous chiral symmetry breaking. 
To perform an extrapolation to the continuum limit we need to account for direct chiral, 
taste (in the staggered case) and rotational symmetry breaking by finite lattice spacing. 

At present, lattice QCD simulations with staggered fermions defined in Ref. [5| reached 
smaller quark masses p compared with other types of lattice fermions. For example, com- 
pared with unquenched simulations using Ginsparg- Wilson fermions (defined in Ref. Q) 
simulations with staggered fermions are computationally less demanding. The computa- 
tional cost depends on many factors like lattice spacing and the practical implementation of 
the fermions in the simulations (see Ref. and references therein). From Ref. |H it follows 
that the simulations with Ginsparg- Wilson fermions may be about ten to one hundred times 
more expensive than with improved staggered fermions at comparable masses. 

However, the price of these advantages of staggered fermions is that they contain an 
internal flaw. In those simulations, to reduce the number of taste degrees of freedom, the 
fourth-root trick is employed. In the continuum limit this trick is consistent because there 
are no taste-changing interactions, however, when one includes these interactions the finite 
lattice spacing makes the root of the determinant non-local and the trick remains contro- 
versial (see, for example, Ref. 0]). In order to include the relatively large discretization 
effects in sta gger ed calculations, staggered chiral perturbation theory has been developed 
(see Refs. [l^, [01 Q, HH). Application of the staggered Lagrangian allowed one to deter- 
mine the pseudoscalar masses and decay constant [ill EH HI- F° r example, recently in 
Ref. M there has been considerable success in determining the B meson decay constant 
using the MILC collaboration unquenched gauge configurations with three flavors of light 
sea quarks. In these calculations both valence and sea light quarks were represented by the 
highly improved (AsqTad) staggered quark action which allows for a much smoother chiral 
extrapolation to physical up and down quarks than has been possible in the past. 

Being computationally very costly, Ginsparg- Wilson fermions, on the other hand, 
do not suffer from the inconsistency of staggered fermions. Moreover, these fermions 
have an exact chiral symmetry. It has been shown in Refs. 0, 0] that it is possible to 
combine the good chiral properties of Ginsparg- Wilson fermions, taking them as valence 
quarks where they are computationally less demanding, with numerically fast staggered 
fermions as sea quarks. Extrapolation for the hybrid (mixed fermion) program has 
recently been addressed for baryons by Tiburzi in Ref. 0| . Here we consider the application 



of this program for the p meson, which is a classic testing ground for QCD with light quarks. 



In this report we extend the study of QCD for staggered sea and valence Ginsparg- Wilson 
quarks in Ref.jil to calculate the mass of the vector meson at one-loop level up to order 
0(a 2 ). QCD with staggered sea and valence Ginsparg- Wilson quarks, which as in 0, 0] 
we will refer to as mixed QCD, consists of N v fully quenched (valence) quarks (with N v 
ghost-quarks) and 4Nf staggered (sea) quarks. The masses of the sea and valence quarks, 
in principle, could be chosen to be completely arbitrary, except that the quenched quark 
masses are fixed to be equal to the masses of their ghost partners. The full graded chiral 
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symmetry of mixed QCD in the chiral and continuum limits is described via the semi-direct 
product Gmqcd = SU(AN f + N V \N V ) L x SU(AN f + N V \N V ) R x U(l) v , where AN f and jV„ 
are the numbers of unquenched (staggered) and quenched quarks, respectively. As in full 
QCD, we expect that the axial U(1)a symmetry is also broken by the anomaly We also 
phenomenologically include in our calculations the one-loop contribution from the p — > titt 
decay which cannot be consistently incorporated in the heavy meson formalism because of 
the high momentum involved in this process. However, the inclusion of this loop is justified 



by the physical results obtained, for example, in Refs. |12j, 1121 Il8[ and references therein. 



SYMANZIK ACTION 

In order to consistently include lattice discretization effects the Symanzik action [l9l . |2 
needs to be constructed, based on the symmetry constraints of the underlying lattice theory, 
with further projection to the low energy effective theory. This task was carried out in 
Refs. (3), 0| and here, we will only briefly mention the main results. 

In general, the Symanzik action up to order 0(e A ) has the form: 

S S ym = S P q + aSi + a 2 S 2 + a 3 S 3 + a 4 S 4 + ... (1) 

where Sp n is the action corresponding to continuum partially quenched theory described in 
Refs. 21, 2^1, S n (n G N) is the action consisting of dimension (n + 4) operators. Here, 



we will require that the heavy meson chiral Lagrangian, corresponding to Spq, matches 
the corresponding partially quenched chiral Lagrangian, similar to that in Ref. [23| or in 
Ref. Il7f when the lattice spacing is set to zero. From discussions in Ref. (see also 



Refs. |2J, |26|, |28[ ) the actions Si and S3 are equal to zero because of the specific character of 
the fermions we use 1 . The action S2, according to Ref. 0, contains only six allowed operators 
which break taste or rotational symmetries (see also for details Refs. @, 0, HHj] ) . The chiral 
Lagrangian corresponding to the action S^ will give only a tree level contribution of order 
0(a 4 ). Although allowed in our counting scheme, this contribution will be disregarded, as 
currently we are interested in the 0(a 2 ) analytic behavior of vector meson masses. 

Using the results from the operator construction in Refs. 0, |ij we will study the ana- 
lytic behavior of the vector meson mass in both the chiral and continuum limits, without 
constructing the explicit expression for the heavy meson PQChPT. 



GOLDSTONE MESON MULTIPLET SECTOR 

Mixed ChPT in the Goldstone meson sector with staggered and Ginsparg- Wilson quarks 
was first studied in Ref. [3]. We first give a brief description of that work. 

Assuming that the Gmqcd symmetry is spontaneously broken down to its vector part, 
the particle spectrum will contain light pseudoscalar bosons. These Goldstone meson fields 
can be written in terms of a (4Afy + 2N V ) x (4Nf + 2N V ), unitary matrix field, S, defined as: 

£ = exp(2z$//), (2) 



There are no dimension 5(7) operators which can be consistently formed by the quark bilinear and which 
preserve the symmetries of the theory in the approximation we are interested in, see Ref. 0| and references 
therein. 
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where £ G U(4Nf + N V \N V ), $ is the bosonic matrix and / is a low energy constant. In 
particular, we are interested in the case with Nf = 3 and N v = 2, for which the bosonic 
matrix has the following form: 
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Following the notation of Ref. 
X —* xx, Y 



quarks), R% x — > xx, similarly for R xy , Ry X , and Ry y 



the following correspondence is implied: P — > xy, 
yy, {P + , X and Y are the analogous combinations of valence ghost 

Qfv — > Fv, where F G {u,d,s} and 
i> G {a;, y} (note that Qf^ is a 4 x 1 matrix in taste). The pseudoscalar states U, D, 
etc. in Eq. (jSJ), which consist of staggered quarks only, form the "sta ggere d" sector of the 
mixed pseudoscalar boson ChPT. This sector is well studied in Refs. jRMlll fl5| . Each state 
in the "staggered" sector can be represented as: 



U 



16 rp 

5>T. 



F a — {^5, ^5; ^fiw, C/i, Clj 



(4) 



a=l 



(similarly for D, K + , etc ) where F„ are the sixteen Euclidean gamma matrices associated 
with taste (see, for example, Ref. |13j). 



Under chiral symmetry transformations: 



(5) 



where L{R) G SU{AN } + N V \N V ) L{R) . As in Refs. @, Q, we adopt the following counting 
scheme: 



m, 



2 a2 rv> 2 
" IV QCD — fc > 



(6) 



where Aqqd denotes the typical QCD scale. This counting scheme is relevant for simulations 
with improved staggered quarks |27( . The leading order chiral Lagrangian (in this counting 
scheme) that contains the terms of order 0(p 2 ,m q ,a 2 ) is of the form: 



r- 



fB 



777 

< SM f + E+M > + — ^ < $ 2 > + a 2 V 



(7) 



4 6 

Here < ... > denotes a supertrace in flavor space and B is a low-energy constant. In our 
case the diagonal mass matrix is: 



M = diag(m u £i, m d £i, m a £j, m x ,m y ,m x ,m y ). 



(8) 



For convenience, to derive the flavor neutral propagators, the singlet mass parameter m 2 , in 
Eq. (JZJ) is allowed to be finite, however, at the end of the calculations we always could take 
the limit ml — > oo. 
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The potential V in the leading order Lagrangian comprises all terms proportional to a 2 . 
At this stage we will not go deep into describing the form of the V (we refer to Refs. 0, 0| 
for more information). Here, we present only those results required for our further study. 

At tree level, because of the specific properties of potential V, the valence-valence mesons 
with only Ginsparg- Wilson quarks obey the continuum-like mass relations: 

m 2 p = B(m x + m y ), m 2 x{y) = 2Bm x ^ y) . (9) 

In our calculations the only relevant valence- valence Goldstone bosons will be flavor neutral 
because the valence quarks in the sea loops will be canceled by their ghost partners and, 
in addition, G-parity suppresses the process 2 p — > tj'tt. We also will be interested only in 
flavor-neutral and taste singlet staggered quark combinations, so in the sea-sea sector, for a 
meson in a singlet taste channel / made of sea quarks of flavor F, we have: 

m 2 Fj = 2Bm F + a 2 A! (10) 

where the constant Aj is given, for example, in Ref. |2| (Eq. (27)). Finally, the masses of 
the valence-sea mesons (Fv or Fv) are given as: 

m Q Fx{y) = B ( m F + m x(y) ) + a 2 A Mix , (11) 

where A Mix is an undetermined parameter in this mixed theory which is taste independent. 

The flavor-charged (non-diagonal) fields in Eq. ((SJ) have only connected propagators 
(G™ n ), in the quark- flow sense, which are given by the expression 3 : 

G ^ - 7fk (12) 

where i, j G {ghost, quark}, £j = 1 if % = quark and — 1 if i = ghost. Disconnected propa- 
gators for flavor neutral mesons (like, for example, U or X) contain hairpin-like interactions 
(throu gh t he Wq term). The sum of all these hairpin-like diagrams (in the way described 
in Ref.|l3l|) for the disconnected XY propagators (G dlsc [XY\) will give (as in Ref. ([3]) the 
following result: 



Qdisc \J£Y] — m ° — ' "~ U I' X ' 1 ' "~L>li^' 1 1 ""Si) (13") 

3 (q 2 + m 2 x )(q 2 + m\)(q 2 + m 2 o)(q 2 + m 2 j)(q 2 + m 2 ^) 

where 7r° , rjj and v( l are the diagonalized flavor-neutral, singlet taste-channel mass eigen- 
states. Although there should be vector and axial-vector hairpin contributions, as in 
Ref. [13], we do not include them here because they are irrelevant in our later calculations. 



2 (q 2 + m 2 Ul ) (q 2 + m\ ) (q 2 + m% ) 



MIXED ONE-LOOP CORRECTIONS TO THE VECTOR MESON MASS 



ChPT with staggered sea and chiral (Ginsparg- Wilson) valence quarks was successfully 
applied to calculate the corrections to pseudoscalar masses and decay constants to one 



2 p and rj' have even G-parity, while the pseudoscalar meson has odd G-parity 

3 we will work only with Euclidean propagators 
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loop in Ref. [3( and recently, in Ref. J4J, this approach was applied to calculate masses 
and magnetic moments of baryons. In this work we will calculate vector meson mass 
corrections at one-loop level, without explicitly constructing the corresponding heavy meson 
Lagrangian but rather using general principles. 

From Eq. it follows that continuum heavy meson chiral Lagrangian will consist of 
the parts: L x ~ 0{p) ~ 0(e) and L 2 ~ 0(m q ) ~ 0{e 2 ). As shown in Ref. [TJ, from the 
comparison with the simulations in Ref. the contribution from the terms of order 0{m 2 q ) 
may be neglected 4 . The chiral Lagrangian, Lb, corresponding to S 2 in Eq. (Q), will consist of 
terms of order 0(a 2 ) ~ 0{t 2 ). And finally, we will not have terms of order 0(p n+1 ) (n e N) 
in the heavy meson chiral Lagrangian because they can be eliminated by using the equations 
of motion [29|. 

Now, it is easy to see that L x gives one-loop chiral corrections to the vector meson mass 
- see the diagrams (a, c, d) shown in Fig. 1. The terms in L 2 and will give tadpole as well 
as tree level contributions. From these considerations for the total tree level correction Otree 
to the vector meson mass (up to order C(e 4 )) it follows that: 

ow(a*) = cia 2 + c 2 m x + c 3 m y + c A m u + c 5 m d + cam s (14) 

where, in general, q = Cj(/i) (i = 1, 6) with fi as the renormalization scale. 

For concreteness we will calculate the mass correction of a charged p meson at rest. (The 
calculation is somewhat more complicated for a neutral p meson, but because of isospin 
symmetry, the mass correction of all p mesons must be the same.) 




f,x— y 



FIG. 1: The quark flow diagrams 



In the case Nf = 3 and N v = 2 we have staggered u, d and s sea quarks with x and y 
valence quenched quarks (corresponding to different flavors). The quark flow diagrams for 



4 Although Ref. considers effective theory with different types of fermions, the results in the continuum 
limit should match. 



6 



this case are presented in Fig. 1, with the propagators for the connected diagrams (Fig. 1 
a, b) given by: 
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(15) 



where p is the initial momentum of the vector meson with mass M (p 2 = —M 2 in Euclidean 
space), q is the momentum of intermediate pseudoscalar meson state and instead of using 
a heavy meson propagator we just used the usual one (which is practically the same) with 
the mass of intermediate vector meson equal to the mass of the initial one. The subindices 
on the propagators remind us that in the self-energy integral we need to take into account 
the vector-vector-pseudoscalar (vvp) and vector-pseudoscalar-pseudoscalar (vpp) nature of 
the meson interactions. For the diagrams in Fig. lc and d in the limit m — > oo, we have: 
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(16) 



Notice, that for case shown in Fig. 16 there are no additional double hairpin diagrams 
because of G-parity. Now, let us take m u = = m (Nf = 2 + 1 picture) as in MILC 
simulations [30]. In this case: 
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(17) 



Let us also require that m x = m y and mu l = mx or equivalently Bm x = Bm + |a 2 A/. The 
first choice will help us to find the correction to the p meson mass because it consists of two 
light quarks. The second choice is natural because in the continuum limit it gives m x = m. 
Taking these two conditions into account we can write the above propagator in the form: 
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Using the method in Ref. [3l( for the total self-energy we have: 
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(19) 



(20) 
(21) 
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we used the notation u 2 = k 2 +m 2 and the functions u 2 v ,.(k 2 , A) are finite range regulators 
( [Til 31, ll^, Hil). For example, the leading non- analytic (LNA) contribution to the self 
energy is: 



« = -f^ Um% ux + 2m% sx - 2m x + |< ) i 22 i 



In the continuum limit (where we have m x = m) the LNA contribution is: 



yLNA n 2 / 9 

sm lna = ^ = -^-(2ml + 2m% + ) . (23) 



This result agrees with that in Ref. |34j for the vector p meson - where in the limit N c — > oo 
they found: 



6M ? = "12^ 



2m^ + 2tt& + jjmj 



(24) 



Note that while the coefficients in these two expressions are written differently, since g = 
2g 2 /f for values of g = g pu)7 , = 16 Gel/" 1 (0), g 2 = 0.75 (0) and / = 93.75 MeV, the 
expressions are equal. 

It would be also interesting to consider the case, when m u = m d , mx = my I and 
my = iris i { m y — m s in the continuum limit). Using these conditions and recalling that 
the K* vector meson consists of one light and one strange quark we find the self energy 
contribution: 

Ztotai = ( £ ™ P ( M >^, A) + S OTp (M 2 ,mJ /s ,A) + S^fm^^m^, A)j - (25) 



f=u,d,s 

-]-Z vvp (M 2 ,m x ,A) - Z vvp (M 2 ,m 2 Y ,A) + ~ 
2 6 



In the continuum limit (m x = m u and m y = m s ) the LNA contribution to the total self- 
energy is: 

^LNA _ 92 /3 3 . Q a . 1 3 



E ^ = i^(i< +3 < + i m V' (26) 

Again, comparing this result with that in Ref. [34| for the K* meson mass correction (in the 
limit N c — > oo): 



-m% + 3m 3 K + -m* 



(27) 



we see good agreement with our result for the LNA contribution in the continuum limit. 

Now, returning to p meson, let us take into account the tadpole contribution from L 2 (&) ~ 
0(a 2 ,m q ). This case is similar to L\ but without intermediate vector meson propagators 
and we don't need to include the diagram with two intermediate pseudoscalars. The total 
self-energy correction from i s: 

% otal = 4Z tad {m 2 Qux , A) + 2Z tad (m 2 Qsx , A) - 2E tad {m x , A) + ^ tad {m 2 m , A), (28) 
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with the definition: 



v, / 2 *n n f dAk k 2 u 2 (k 2 ,A) 

where, in general, 3 = J3'm 2 + 3a 2 , 3' and /3 are undetermined parameters and u(k 2 ,A) is a 
finite range regulator 0, I32I l33f . The term proportional to 3' is absent in the continuum 
limit (see, for example, Ref. |l7( when lattice spacing is set to zero). This means that the 
tadpole contribution comes only from the operators in L&. Then, for example, if u — 1 
(applying dimensional regular ization) we have: 

S iad (m 2 ) = ^- 2 m 2 ln^ ~ 0(e 4 ). (30) 

The calculations above were performed in the case where x and y correspond to different 
flavors. To calculate the mass of the neutral vector meson we need to change our quark flow 
diagrams and add one more with DHP insertions on each intermediate pseudoscalar line. 
However, as noted above, isospin symmetry for the vector meson, means that we need not 
perform separate calculations for the neutral vector meson mass. 



VECTOR MESON MASSES UP TO ORDER 0(e 4 ) 

Imposing the conditions m x = m y = m u + ^a 2 A/ with m u = m<i we can simplify the 
expression for the tree level contribution and the vector meson mass in both chiral and 
continuum expansions can be written in the form: 



My 



(c + c a a 2 + c q m 2 n + c s m 2 S[ ) 2 + T, total (/i) + S tota/ (/i), 



(31) 



where m 2 = 2Bm x = m 2 x , m 2 s = 2Bm s +a 2 Ai and c, c aj(?iS are new dimensionful parameters. 
The dependence on m 2 Si may be absorbed in the coefficients c and c a if one fixes the mass 
of the strange quark in the simulations. 

From the definition of T, vvp and T, vpp using the appendix in Ref. 3] we can separate the 
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where u p is the physical mass of the p meson and using Eqs. (JOHHJ) we can express all meson 
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masses as functions of m\ and m\\ 



m\ ux =ml + a 2 \A Mix + -A, J , (33) 

,2 1 2 , -^ 2 i „2 ' • ' 



m m = l m * + \ m2 Si- ( 35 ) 

It is also of practical interest to consider the general case with only two conditions: m u = m^ 
and m x = m y . In this case the total self-energies are: 



f=u,d,s 

4 4 
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a= l 



a 



{m 2 Vi - m 2 x ) 



1 - 



m\ - m\ ' 



from which one can see if m x = m n = mjjj then we have the same result as in Eqs. (|25J28|) . 

To perform chiral-continuum extrapolation of lattice data for the p meson with Ginsparg- 
Wilson quarks in a staggered sea the consistent numerical evaluation of self-energy integrals 
^totaiin) and Etotaz (/-0 nee ds to be done using the particular type of finite range regulator. 
The extrapolation formula, in a numerical simulation with fixed strange quark mass, will 
have five coefficients, some of which are functions of the FRR cutoff. 



CONCLUSION 



Based on the recently developed mixed PQChPT for pseudoscalar meson sector we have 
extracted a chiral-continuum extrapolation formula for the vector meson masses at one-loop 
up to order 0(a 2 ). Instead of writing explicit expression for the vector meson Lagrangian 
we have used general principles in order to calculate the mass corrections at one-loop order. 

Up to order C(e 4 ), vector meson masses, in general, depend on six, presently unknown, 
parameters: c, c a>q>S) (3 and A Mix . The low energy coefficient A Mix is model independent and 
may also be estimated from the extrapolation of pseudoscalar masses. If one, however, is 
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interested in vector meson mass corrections up to order 0(e 3 ) then the number of parameters 
reduces by two: 

Ml = {c + c a a 2 + c q ml + c s m 2 Si ) 2 - M- (2m 3 Qtm + m 3 Qsx - m\ + X -m\\ . (36) 

This can be used to test the theory predictions up to order 0(e 3 ). It is also possible to 
reduce the number of extrapolation parameters by one if we fix the strange quark mass in 
the simulations. 

During the derivations the following approximations were used: m x = m y , m u = and 
mx = vnu I for the p meson. This approximation may be adjusted, for example, to the K* 
vector meson if we assume m u = m^, mx = Tnjjj and my = rrisj- The idea beyond this is 
that in the continuum limit the masses of the valence quarks, m x and m y , are required to 
correspond to the masses of valence quarks from which the vector meson is constructed. 

In our calculations we include (phenomenologically) the relativistic p — ► mi decay 
process, which although cannot be consistently included in the framework of low energy 
effective field theory, is required to give consistent physical predictions. Finally, we did not 
take into account the finite volume effects from the lattice, which are already well studied 
in Refs. |I1 B H . 
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